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Abstract 

All the orbital Ml excitations, at both low and high energies, obtained from 
a rotationally invariant QRPA, represent the fragmented scissors mode. The 
high-energy Ml strength is almost purely orbital and resides in the region of 
the isovector giant quadrupole resonance. In heavy deformed nuclei the high- 
energy scissors mode is strongly fragmented between 17 and 25 MeV (with 
uncertainties arising from the poor knowledge of the isovector potential) . The 
coherent scissors motion is hindered by the fragmentation and B(M1) < 
0.25 fj,pj for single transitions in this region. The (e, e') cross sections for 
excitations above 17 MeV are one order of magnitude larger for E2 than for 

Ml excitations even at backward angles. 
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I. INTRODUCTION 



The nature of the orbital magnetic dipole (Ml) strength in deformed nuclei is still a 
controversial matter. The low-lying orbital excitations with = 1 + have been interpreted 
[ lj as isovector rotational vibrations, i. e. as a manifestation of a weakly collective scissors 
mode ||. A strongly collective scissors state was predicted by the classical picture of the 
Two-Rotor Model ||, assuming that neutrons and protons perform out-of-phase rotational 
oscillations around an axis perpendicular to the nuclear symmetry axis. The fragmentation 
of the Ml strength over a broad energy region, found experimentally between 2 and 4 MeV 
in heavy deformed nuclei (e. g. Refs. fU) and present in microscopic theoretical calculations, 
hinders the comparison with the simple classical picture. It has been shown earlier |l]g 
that the energy and the amount of the high-lying Ml strength, corresponding to AJ\f = 2 
excitations, are strongly dependent on the type of the residual quadrupole-quadrupole (or 
non-spin) interaction used in microscopic calculations within the quasiparticle random-phase 
approximation (QRPA). 

The above interpretation in terms of scissors mode is supported by some theoretical 
works [|6|-|T0|1 and questioned by others [fLl-151- The scissors mode is found in Refs. [|i~3| -|T5 



at high energy, E theor « 22 — 24 MeV, as a part of the isovector giant quadrupole resonance 
(IVGQR). These three works use Landau-Migdal interactions in QRPA. Similar conclusions 
are obtained in QRPA with separable forces [|l(J, where the IVGQR lies at E t h e or ~ 27 MeV 
and overlaps 40% with the scissors mode. It was shown [14[] that the collective hydrodynam- 
ical model predicts a high-energy scissor mode, but the low-lying state is spurious. However, 
after taking additionally the nuclear elasticity within the same classical model into account, 
a low-energy mode was obtained as well ||15|| . The isovector rotor model |l7j predicts the 
appearance of the scissors mode at both low and high energies. 

We are going to discuss in Sec. Ill the isovector coupling of this interaction, which 
influences strongly the high-lying orbital Ml strength. An isoscalar symmetry restoring 
interaction is introduced in Sec. II, together with an isovector interaction. The E2 strength 
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distribution is studied in Sec. Ill in relation to the isovector coupling constant. The high- 
energy Ml strength and related (e, e') cross sections are discussed in Sec. IV. The conclusions 
are summarized in Sec. V. 



II. DECOUPLED ISOVECTOR INTERACTION 

The deformed mean field Ho is obtained in our formalism by diagonalizing an 

axially-symmetric Woods-Saxon potential. Electric and magnetic transition operators are 
treated on equal footing by introducing a signature m — ±1 0, corresponding to the 
directions x and y, indistinguishable in the intrinsic frame. For instance, the total angular 
momentum operator, J t = J(m,t), defined by Eq. (4) of ||, has quasiparticle (q.p.) matrix 
elements j l ki = j(ki,mt) The isospin index t denotes neutrons or protons. We shorten 
the notations in the following by omitting, wherever possible, the signature index m, which 
is of technical interest only. 

Following the procedure of Pyatov PU| , we construct in the quasiboson approximation 



operators of quadrupole type F. They are used to define the isoscalar interaction H5 
which restores the rotational symmetry of the mean field H , violated by the deformation: 



F t = F(m,t) = [H ,J t ], 
Fs = F n + F p , F v = F n — F p , 

2 

Under the assumption for a simple quadrupole deformation Q 20 only, as in the case of the 



H s = -^Y>F|(m), (1) 



Nilsson potential, the operators F s (m) ([[]) are proportional to the quadrupole operators 
Q2.ii 

A separable residual interaction of the general form, 



H ^ = -\ E m [ fc » F «( m ) + k PP F l 

m 

(m)Fp(m) , (2) 
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can be decomposed into a sum of isoscalar, isovector, and coupling terms with the corre- 
sponding coupling constants, 

k\ ~^(k nn -\- kpp 2/c n p), 
k = ^(k nn — kpp), r = ki/k$. (3) 

The isovector operator Fy ([j]) is obtained from the isovector angular momentum J n — J p . 
However, it is more convenient to construct a relative isovector interaction Hr and the 
coupling He with the isoscalar term (0) by introducing a relative angular momentum J R : 



^■Jn-y^Jp' X t — F t ]), 



F R = F R (m) = [n ,J R {m)), 

H c = -^rE™ F *(™) F *M- (4) 



It is easy to show that the relative isovector interaction Hr is rotationally invariant and its 
coupling constant k R remains a free parameter. This interaction is slightly more symmetric 



than that used in Ref. ||12|| , but the difference is insignificant. The condition for rotational 
symmetry provides the same value for the isoscalar constant ks as in Eq. (9) of || and 
closes the coupling He between the isoscalar and relative channels: 

[H + Hs + Hr + H c , J] = 0, 

=► k s = 1/X, k c = 0, X = X n + X p . (5) 

Thus, the use of relative interaction enables the decoupling of the two channels and the 
symmetry-restoring procedure is no more necessary, since the hamiltonian (|) is already 
rotationally invariant. 

We introduce further a spin-spin residual interaction H55 ||18|| , which is always rotation- 



ally invariant and does not change the coupling constants (|5]). The latter can be related to 



the coupling strengths of the quadrupole interaction in the general form (0). Thus, our q.p. 
model hamiltonian, used in the present calculations, has the following form: 



H = H + Hff + H55, 

X X 
Kn = k s (l + r—r-), k pp = k s (l + r—^-), 

k 

k np = k s {l-r), r = —, (6) 

where k$ is calculated microscopically (|5|) and the ratio r is treated as a free parameter to 
be discussed in the next section. The relationships with the constants k , fci, and k! are 
obtained upon insertion of (||) into (|3|). The constants of the spin-spin interaction, 

c(+)A = 200 MeV, c(-) = -0.5c(+), (7) 



are derived from nuclear matter calculations in the way described in Ref. ||21|| . The 
deformation-independent Woods-Saxon parameters are taken from Ref. ||22||. When applied 



to other rare-earth |23| and actinide [^] nuclei, this parametrization has lead to a good 
agreement with experimental data for single Ml transitions at low energy, observed in (e, e') 
and (7, 7') experiments. The energy distribution of the spin-flip Ml strength between 6 and 



10 MeV, determined by inelastic proton scattering, has also been reproduced in f21]j23 |. 

The QRPA equations of motion are solved for if" = 1 + excitations in the intrinsic frame 
using the model hamiltonian (|j). We present here as an example the results for 160 Gd, 
obtained with deformations (3% = 0.26, (3& = 0.025 and pairing gaps A n = 0.8 MeV, A p = 
1.1 MeV. Expressions for Ml and E2 transition probabilities in terms of RPA amplitudes 
are given in Ref. [0]. The B(M1) values for transitions to I n K = 1 + 1 states are calculated 
with bare orbital and effective spin gyromagnetic factors, g s = 0.7gf ec . The B(E2) values 
for transitions to the corresponding I n K = 2 + l states are obtained without introducing 
effective charges, i. e. they are of purely proton nature. 

Our attention shall be focused in the next section on the ratio r determining the 
strength of the isovector spin-independent interaction. 
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III. ISOVECTOR COUPLING CONSTANT 



The energy distribution of the E2 strength with K = 1 is shown in Fig. 1 as histograms 
for three different values of the isovector coupling constant k^ (||). The latter is fixed by 
the ratio r = k^/ks, because the isoscalar constant k$ is determined microscopically from 
(H). The lower bump has a peak at about 11 MeV, which is smeared out by the histogram. 
It represents the K = 1 component of the isoscalar giant quadrupole resonance (ISGQR), 



observed by inelastic electron scattering [24]. Its position, determined by the isoscalar 



constant k$, agrees with the prediction of the vibrating potential model |25| 



E(21;is)^^ I - 3 (l- 5 -)MeV. (8) 
The summed contribution of the K = ±1 components to the classical isoscalar and isovector 



energy-weighted (e.w.) sums are |26 



S(E21-is) c = S(E21) c Z/A, 
S(E21; iv) c = S{E21) C N/A, 
S(E21; is) c + S{E21; iv) c = S(E21) C , 

We calculate microscopically the mean-square-radius, (r 2 ) p = 27.3 fm 2 , and obtain with 
this value from (|9|) the classical e.w. sum S(E21)c = 57573 e 2 /m 4 MeV. The RPA 
isoscalar E2 e.w. strength between 9 and 11 MeV, shown in the middle plot of Fig. [I], 
is S(E21; 9 — 11 MeV) t h = 12215 e 2 /m 4 MeV, i.e. only 53% from the classical isoscalar e.w. 
sum S(E21;is)c = 23029 e 2 /m 4 MeV are exhausted. This result corroborates with the ex- 
perimental data for other heavy deformed nuclei ||24|| . However, no strength is missing, since 
the total RPA e.w. strength, S(E21; 0-30 MeV) th = 57272 e 2 /m 4 MeV, exhausts 99.5% of 
the above classical value S(E21) C - One should note further that the decomposition (Q) of 
the classical E2 strength into isoscalar and isovector parts is a rather crude approximation, 
because it does not take into account the redistribution of these two components caused by 
the residual interactions. 



It is seen from Fig. [I] that the low energy E2 strength distribution and the position 
of the ISGQR are rather insensitive to the isovector residual interaction. The latter is 
repulsive and concentrates the isovector strength in the higher-lying bump, consisting of 
AjV = 2 excitations. The non e.w. E2 strength at high energy decreases for a stronger 
isovector interaction, because the e.w. E2 strength is conserved in the three plots of Fig. 
|l|. This is due to the fact that the interaction (§) redistributes the E2 strength, but does 
not contribute to the E2 e.w. sum rule. The latter is determined by the mean field and the 
spin-spin interaction which are fixed in the three cases of Fig. [l|. 

The isovector coupling constant could be determined from the position of the second 
maximum, representing the isovector giant quadrupole resonance (IVGQR). Experimental 
data on it for heavy deformed nuclei are scarce ||27|| . The measured (e, e') spectrum is usually 



fitted by a number of gaussians. That one, lying closest to the theoretically expected energy 
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E(21;iv)*^(l- 5 -)MeV, (10) 

is identified with the IVGQR. The prediction ( |10"D results from evaluation of the isovector 
coupling constant within the collective model. The main uncertainty of this procedure 
originates from the isovector potential V±, which is estimated by relating the isovector density 
with the potential asymmetry energy 

P A 

Vx = 8K as . (12) 



Using a recent empirical value, Ky = 30.6 MeV |28|], for the volume asymmetry energy from 
the Bethe-Weizsacker formula, and subtracting the kinetic energy bki n = 2Kki n = 25 MeV 
one obtains K as = Ky — K^in ~ 18 MeV. The resulting isovector potential V\ = 144 



MeV ( 12]) is a little bit larger than the previous estimate of 130 MeV |2{J , based on a smaller 
value of Ky. 

The identity ( |T2"D follows from ( |Tl~D under the assumption that the integral in the l.h.s. 
of (|TT|) is equal to (N — Z) 2 /A. This is true for constant neutron and proton densities 



within a sphere with a sharp surface. When the integral is calculated with more realistic 
Woods-Saxon densities and the isovector potential is determined from ([11]), instead of (|T2"|). 
we obtain V\ = 67.5 MeV for 160 Gd. Test calculations for several more rare-earth nuclei 
have lead to values close to 68 MeV. 

The ratio between the isoscalar and isovector coupling constants of the quadrupole in- 
teraction can be calculated within the collective model p5| . Using our microscopic mean 



square radius for 160 Gd, the above isovector potential V\ = 67.5 MeV, and an oscillator 
constant TiloqA = 41 MeV, we obtain for this ratio r' = k\/ko = —1.5. However, as discussed 
in Sec. II, the operators F (j|D do not coincide with the quadrupole operators. Moreover, 
the isoscalar coupling constant k$ is already determined from the rotational invariance (^j). 



The isovector constant k R can be evaluated in our case following the procedure of Ref. |25 
for the operator F: 



f F f pFdr = -(FF" 1 ") -1 , 

J J 8 

(FF f ) = 2j2 E ikJik> (13) 

ik 

where the summation runs over neutrons and protons. An additional factor of 1/2 is present 
in flT3| ) in comparison with [^J , since our coupling constants (0J^) are two times smaller than 
usual. This is due to the m-symmetrization of the interaction Hff (§,§)• Using V\ = 67.5 
MeV, one obtains from (|5|JT3|) the ratio r = k^/ks = —1.67, consistent with the above value 
of r' for a quadrupole interaction. 

These numbers contain still rather large uncertainties. Corrections to the isovector po- 
tential could arise, e. g. from the surface asymmetry energy. It contributes even to the 
IVGQR regarded as a predominantly volume oscillation. Further renormalization ef- 



fects []3TJ can lead even to attractive n-n and p-p quadrupole interactions with r' = —0.6. 
Thus, the theoretical estimates of the isovector coupling constant are not accurate enough 
and only the experimental data on the IVGQR could decide on this issue. We adopt the 
value r = —2, producing with our potential an IVGQR centered at about 22-23 MeV, as 
seen from the middle plot of Fig. |l]. 
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This result agrees with experimental findings in several rare-earth nuclei |27] and more 
reliable data on 208 Pb ||24j| . The theoretical E2 strength between 17 and 25 MeV, shown in 
the same plot and corresponding to the IVGQR with K* = 1 + , amounts to 1475 e 2 fm 4 . 
The value is close to a rough estimate [0 from a schematic two-level model and agrees 
qualitatively with recent results on the IVGQR in 154 Sm, obtained from inelastic scattering 



of polarized protons f32[. The microscopic isoscalar coupling constant ks (fj) gives rise to 
the isoscalar GQR with K n = 1 + at 11 MeV. Its E2 strength from the middle plot of Fig. 
H, summed between 9 and 11 MeV, is B(E2) = 1193 e 2 /m 4 , 

in agreement with experiment 

m 



IV. HIGH-ENERGY Ml STRENGTH AND (e, e') CROSS SECTIONS 

The energy distribution of Ml strength in 160 Gd, summed in bins of 1 MeV, is shown 
as histograms in Figs. |2] and [3] for different isovector ratios r (|6[). The total and orbital Ml 
strengths are represented by contour lines and shaded areas, respectively. It is seen that 
the high-energy Ml strength is purely orbital, while the low-energy one contains spin con- 
tributions, which are dominant between 6 and 10 MeV. Comparison with the corresponding 
cases from Fig. [l| shows that the high-lying orbital Ml strength resides in the same energy 
region as the E2 IVGQR. 

The low-lying orbital Ml strength is influenced only slightly when increasing the mag- 
nitude of the isovector interaction from r = to r ~ — 1. A further increase (r < —1) does 
not produce appreciable changes. The quadrupole interaction (H) contributes to the Ml 
e.w. sum-rule Jl[. The repulsive isovector interaction shifts the high-lying orbital strength 
to higher energy and increases the e.w. sum. 

Figure displays results corresponding to the adopted isovector coupling constant, r = 
—2. The lower plot of Fig. ^| was obtained by including only 2q.p. configurations below 20 
MeV in the RPA calculations, as done also in [TJ|. A schematic, single excitation is produced 
at 22 MeV, carrying most of the high-energy strength: B(M1) = 3.6 fij^, B(E2) = 
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1305 e 2 fm A . However, this strongly collective state is no more present in more realistic 
calculations, including further 2q.p. configurations up to 30 MeV, as seen from the upper 
plot of Fig. |2|. The single excitation is strongly fragmented over a broad energy region in 
this case. 

Comparison with our previous results from Figs. 3 and 4 of Ref. |5[] demonstrates that 
the strong fragmentation is not solely due to the large basis, but results also to a great 
extent from the restoration of rotational invariance. It is seen from these two figures that 
the quadrupole interaction alone produces a high-energy orbital Ml strength, which is still 
strongly concentrated in a narrow peak. In contrast, considerable fragmentation takes place 
after adding the symmetry-restoring interaction [jlj (cases denoted by a prime in Figs. 3,4 of 
Ref. 0). This strong fragmentation is in a qualitative agreement with the present results, 
obtained from the rotationally invariant hamiltonian (|6]). 

The summed Ml strength above 20 MeV in the upper plot of Fig. Q is equal to that of 
the single excitation at 22 MeV from the bottom plot, resulting from the basis cut-off. The 
equality is due to the fact that the high-energy strength originates from 2q.p. configurations 
with AAf = 2, lying at about 15 MeV. Thus, configurations above 20 MeV do not generate 
but only redistribute the high-energy strength, pushed by the repulsive isovector interaction 
from 15 to 22 MeV. 

Let us consider the more realistic results in the upper plot of Fig. |2|. The summed 
(predominantly orbital) Ml strength between 17 and 25 MeV is 3.6 [x 2 N , in agreement with 
a recent estimate for 154 Sm from (p,p') experiments |32]. It is derived from a sum-rule |33 



connecting the Ml and E2 strengths. The strongest single 1 + excitation from the considered 
energy region lies at 22.5 MeV and carries only 0.25 fi 2 N strength, i. e. it could be hardly 
resolved experimentally. This excitation overlaps only 4% with the microscopic scissors state 
|||2|], obtained by rotating the neutron and proton configurations of the deformed mean field 
in opposite directions. The low-lying 1 + excitations below 9 MeV overlap all together 77% 
with the scissors state, while the overlap of those between 17 and 25 MeV amounts to 55%. 
The total overlap of all the RPA excitations below 30 MeV is 1.46. Thus, the low- and 
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high-lying excitations exhaust 53% and 38% of the total overlap with the scissors state, 
respectively. 

Let us note that the total overlap is exactly 100% for a vanishig isovector interaction 
|2||J- The isoscalar channel cancels the spurious (isoscalar) admixtures, but it does not 
contribute to the (purely isovector) scissors mode. Thus, an overlap of 100% is produced 
by the mean field alone. The overlap is not conserved in RPA since it corresponds to a non 
e.w. sum, while only energy-weighted sums are conserved in RPA. The excess beyond 100% 
is generated by the isovector residual interaction, which is not present in the microscopic 
scissors state, constructed from the mean field alone. The overlap reaches 174% in the case 
r = —3.6, plotted in Fig. |[ 

A more meaningful geometric interpretation of the orbital 1 + excitations is provided by 



the comparison with low collective synthetic states [33]. A synthetic state is constructed for 
each RPA excitation by rotating out-of-phase neutron and proton configurations weighted 
with the RPA amplitudes of this excitation. A large overlap of a given RPA state with its 
synthetic counterpart is an indication for an isovector rotational motion, irrespectively of 
the low collectivity caused by the fragmentation. The strongest (low-energy orbital) Ml 
excitation overlaps usually 80-90% with its synthetic counterpart P|,pTJ] . In contrast, no 
one of the strongest orbital states at high energy (around 22 MeV in the upper plot of Fig. 
ID) overlaps more than 7% with its synthetic counterpart. The coherence is restored in the 
schematic single excitation at 22 MeV in the lower plot: it overlaps almost 90% with its 
synthetic counterpart, confirming the presence of the scissors mode at high-energy. 

It is clear from the above considerations that the scissors mode fragments over both the 
low- and high-energy orbital Ml excitations. This holds exactly for the isovector rotational 



model [T7|] , formulated in the schematic two-level basis of the quantized deformed oscillator. 
It was found within this model that the scissors mode exhausts the whole non-spurious 
orbital Ml strength at low and high energy, corresponding to excitations with AJ\f = 0,2, 
respectively. 

A strong isovector quadrupole interaction, based on the operators Q2,±i, leads to similar 
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results [[|, apart from concentrating in a narrower energy region much more high-lying Ml 
strength. Most of it turns out to be spurious. After restoration of the rotational invariance, 
the high- lying strength is reduced and fragmented over a broad region [0,0, in agreement 
with the present results, based on the rotationally invariant hamiltonian (|6]). 

Results from different QRPA calculations are compared in Table I. Our results (first row) 
are similar to those obtained with a zero-range Landau-Migdal interaction |T3l]l4j1 (second 
row), apart from a lesser fragmentation in the latter case. This is seen on the example 
of 164 Dy 0], where the strongest high-lying 1 + state (at 21 MeV) has a large transition 
probability, B(M1) = 1.26 fi%. It seems that the concentration of strength is caused by the 
reduced basis of 2q.p. configurations adopted in Ref. |L4|]. As discussed above, the total 



strength is not influenced substantially by this procedure. 

The results listed in the last row of Table I |16j are similar to those in the bottom plot of 



Fig. |3], because the same basis cut-off (20 MeV) and almost the same isovector ratio are used 
in both cases. The B(M1) value of the single excitation at 26 MeV in the bottom plot of Fig. 
[3] is four times larger than the value 1.5 /i^r, listed in the last row of Table I. The difference 
is due to the fact that 152 Dy, considered in [jnj, has a very small deformation, producing 
less orbital Ml strength than the well-deformed 160 Gd, studied in the present work. 

It is seen from the middle plot of Fig. ||| that a considerable fragmentation takes place after 
including all the 2q.p. configurations up to 35 MeV in calculations. The strong collective 
excitation disappears and its strength is distributed over a wide energy range. As discussed 
in Sect. Ill, the present status of experiment on heavy deformed nuclei does not allow a 
reliable discrimination between the results in Figs. and [3], obtained with different isovector 
coupling constants. 

Theoretical (e, e') cross sections , calculated in DWBA |35| using our RPA wave functions, 
are plotted in Figs. f| and |5]. E2 (dotted curves) and Ml (dashed curves) contributions to the 
cross sections for a scattering angle 6 = 165° are plotted versus the electron incident energy, 
together with their sum (continuous curves). The (e, e') cross sections of the two strongest 
high-energy Ml transitions are plotted in Fig. |4j. Case A corresponds to the RPA excitation 
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with E x = 21.55 MeV and B(M1) = 0.25 while case C refers to the RPA excitation 
with E x = 20.54 MeV and B(M1) = 0.20 fi 2 N . Although both states are composed mainly 
from configurations with AJ\f = 2, the leading 2q.p. components in the corresponding RPA 
wave functions are neutrons with AjV = 4 in the former state and protons with AjV = in 
the latter state. Thus, the former excitation has a more pronounced volume character than 
the latter one. 

In addition to the above two strongest excitations at high energy we consider also the 
two next strongest high-energy orbital Ml excitations from the upper plot of Fig. |2|. These 
four excitations lie between 20.5 and 23 MeV, have a summed B(M1) value 0.73 fi 2 N , and 
a summed overlap with the scissors state 11%. The sum of their cross sections is displayed 
in the upper plot (case C) of Fig. |^. The E2 contribution to the cross section is dominant 
over the Ml cross section for incident energies below 100 MeV. The diffraction structure is 
smeared out, because the considered four excitations have very different transition densities. 

The lower plot (case D) in Fig. [5] corresponds to the schematic collective excitation at 22 
MeV from the lower plot of Fig. |2|. It is seen from Fig. |^ that even in backward scattering the 
cross section is one order of magnitude larger for E2 than for Ml excitations. The displayed 
cross sections are calculated for 9 = 165°, an angle often used in experiments, but we have 
checked that the same E2 dominance is present even at a fully backward angle. This result is 
particularly interesting in view of the well-known expectation J36| for a strong enhancement 
of Ml over E2 excitations in backward scattering. Such a qualitative expectation is based 
on approximations (negligible excitation energy compared to the incident electron energy, 
PWBA), which do not hold for the DWBA cross section of the high-energy excitation, 
dispayed in Fig. [5|. 

The experimental resolution of single Ml excitations is hindered furthermore by the 
strong fragmentation at high energy. One should expect, therefore, that above 17 MeV the 
experimental cross section of inelastically scattered electrons will originate almost exclusively 
from the E2 IVGQR and not from Ml transitions. In contrast, at low energy the Ml cross 



section dominates over the E2 cross section 37 , at least for low transferred momenta. 
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V. CONCLUSIONS 



The high-lying Ml strength is almost purely orbital. It resides in the same energy region 
as the isovector E2 strength (IVGQR). Experimental data on the IVGQR in heavy deformed 
nuclei are scarce. The interpretation of the (e, e') spectra rely on a rough estimate of the 
expected excitation energy of this resonance, based on the collective model. The main 
uncertainty arises from the poorly known isovector potential, which determines also the 
magnitude of the isovector coupling in the residual spin-independent interaction, used in 
microscopic calculations. We have adopted here a ratio r = — 2 between the isovector and 
isoscalar coupling constants. The isoscalar constant is determined microscopically from the 
condition of rotational invariance. On the example of 160 Gd, considered here, the isoscalar 
coupling gives rise to an isoscalar GQR at 11 MeV with B(E2) = 1193 e 2 fm 4 , in agreement 
with experiment. 

The adopted isovector ratio, r = —2, produces in RPA an IVGQR whose K w = 1 + - 
component is centered at about 22 MeV, but strongly fragmented over the interval 17-25 
MeV. Thus, even the strongest Ml transition from this region has a relatively small B(M1) 
value, 0.25 in comparison with the strong low-lying orbital excitations. The high-lying 
strengths, summed in the above interval, are B(E2) = 1475 e 2 /m 4 , B(M1) = 3.6 fj, 2 N . The 
total RPA energy- weighted E2 strength exhausts exactly the classical sum rule for K — ±1. 

The scissors mode is found to fragment over both the low- and high-energy orbital Ml 
excitations. Although the latter are more collective, the coherence of the scissors motion 
is destroyed at high energy by the increased fragmentation, due not only to the high level 
density, but also to the rotational invariance. Most of the strongly fragmented high-lying 
orbital excitations (17-25 MeV) are not performing well-pronounced rotational vibrations. 
This is seen from their small overlap with synthetic states, where such a geometrical motion 
is enforced by construction. Nevertheless, the orbital excitations from the above high-energy 
region have all together a large overlap, 55%, with the collective scissors state constructed 
from the deformed mean field. The 1 + excitations below 9 MeV overlap 77% with the 
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scissors state. The total scissors overlap of all the RPA excitations below 30 MeV is 1.46. 
The excess above unity is produced by the isovector coupling, corresponding to a stronger 
attractive neutron-proton interaction, which enhances the scissors motion. 

The recent canonical quantization of the isovector rotor in relative conjugate variables fT7| 
provides simple analytical results in the deformed oscillator basis. They strongly support the 
interpretation of all the orbital Ml strength, at both low and high energies, as manifestation 
of the fragmented collective scissors mode. 

The (e, e') cross sections of 1 + states above 17 MeV originate mainly from E2 excitations 
even for backward scattering, while the contributions of the high-lying orbital Ml excitations 
to these cross sections are typically one order of magnitude smaller. 

Thanks are due to Jochen Heisenberg for providing us with his DWBA code. This work 
is supported by the Deutsche Forschungsgemeinschaft (DFG). 
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TABLES 

TABLE I. Comparison of different theoretical QRPA predictions on the summed high-energy 
Ml(j) strength and overlap with the scissors state R^|) || in rare-earth nuclei. 



Nucleus 


E 


E B(M1) 


E EB(M1) 


EI(R-i + )l 2 




[MeV] 


K] 


[MeV (1%) 


[%] 


160Q d a 


17-25 


3.6 


76 


55 b 


164 Dy c 


21-23 


3.6 


79 


30 


152 Dy d 


27 


1.5 


40 


40 



a This work. 

b 38% of the total overlap 1.46. 
c Ref. jn . 

d Ref. nrai. 
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FIGURES 

FIG. 1. Energy distribution of B(E2;0^ S — ► 2 + l) values, summed in bins of 1 MeV, 
for 160 Gd, corresponding to three different values of the isovector strength constant (|6|): 
r = -0.6, -2.0, -3.6. 

FIG. 2. Energy distribution of -B(M1;0+ S — > 1 + 1) values summed in bins of 1 MeV (contour 
histograms with a bright shading) for 160 Gd, obtained with isovector strength ratio (|6|) r = —2. 
Dark-shaded areas: only orbital Ml strength. Bottom plot: only 2q.p. configurations below 20 
MeV are included in calculations. 

FIG. 3. The same as in Fig. but for isovector strength ratios r = —0.6, —3.6. Basis cut-off 
in the bottom plot as in the lower plot of Fig. ||. Note the compressed ordinate scale in the bottom 
plot. 

FIG. 4. DWBA (e, e') cross sections of the strongest (case A) and the second strongest (case 
B) RPA excitations around 22 MeV from the upper plot of Fig. g with B{Ml) values 0.25 and 0.20 
[i 2 N , respectively. Scattering angle 6 = 165°. Contributions from E2(0+ s — > 2 + l) (dotted curves) 
and M1(0+ S — > 1 + 1) (dashed curves) excitations of the same rotational band with K n = 1 + , and 
their sum (continuous curves) are plotted versus the electron incident energy. 

FIG. 5. The same as in Fig. ||, but for the summed cross sections of the four strongest Ml 
excitations around 22 MeV from the top plot of Fig. [2] (case C) and the single state at 22 MeV 
from the bottom plot of Fig. |2| (case D). 
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